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SUMMARY 

The solutions to the problem of the near earth satellite with- 
out drag obtained by applying the von Zeipel method and the 
modified Hansen method are compared. Formulas are derived 
for osculating elements when the modified Hansen theory is ex- 
pressed in terms of orbital true longitude. Differences in the 
arbitrary constants are tabulated. Transformations that relate 
the time element of the two theories are also given. 
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INTRODUCTION 

Widely different theories are often used in computing orbits of artificial satellites. It is of in- 
terest to examine the results of different theories when they are applied to the basic problem of the 
near earth satellite without drag. Of special importance are the major theories of celestial mechan- 
ics introduced by Brouwer 1 and by Musen 2 ’ 3 in solving this problem. 

Brouwer 1 applied the method of von Zeipel to the near earth satellite problem and obtained ana- 
lytic representations for the osculating Delaunay and Keplerian elements. The results are given by 
Brouwer to order j 2 in the elements and j 2 2 in the mean motions, where J 2 is the coefficient of the 
second zonal harmonic of the earth’s potential, and equals the quantity 2k 2 appearing in the articles 
of Brouwer and Musen. Musen 2 - 3 , on the other hand, first modified Hansen's method, then by ap- 
plying it to the same problem of the near earth satellite without drag, showed how to obtain the posi- 
tion of the satellite in a semi-analytic manner to any prescribed order of J 2 . The solution of the 
satellite problem in terms of orbital true longitude by Musen 3 is considered below. 

The results obtained by Brouwer are given in a form convenient for comparison with the results 
of many authors. Indeed Kozai 4 , Garfinkel 5 and others have readily compared their solutions with 
Brouwer’s. However, since Musen’s formulations of the problem are intended to provide numerical 
results of high precision for the position of a satellite, explicit analytic formulations of the pertur- 
bations of the elements do not appear in his articles. For that reason, formulas are given here for 
elements derived from the modified Hansen theory in terms of orbital true longitude so that the re- 
sults of Musen and Brouwer can be compared. 

As would be expected, the differences between the two theories are exhibited in the respective 
choices of the arbitrary constants and in the arguments of the trigonometric terms. The constants 
of both theories are discussed and presented in tabular form. The transformations of the variables 
of the angular arguments are presented. Therefore, when the solutions to the satellite problem are 
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carried out to the same order in J 2 by the methods of Brouwer and Musen, full correspondence can 
be obtained by taking into account the differences in the constants and the angular variables. 

THE OSCULATING ELEMENTS 

The definitions of the osculating elements appearing in Brouwer's article may be found in any 
text on celestial mechanics— for example, Brouwer and Clemence 6 . It is, moreover, a relatively 
simple matter to find expressions for the osculating elements of the modified Hansen theory as ex- 
pressed in terms of orbital true longitude; these differ from the corresponding formulas of the 
modified Hansen theory as expressed in terms of eccentric anomaly given by Bailie and Bryant 7 , 
since the W functions differ slightly. To indicate how such representations of osculating elements 
are derived, we shall now review briefly some of the concepts and definitions of the modified Hansen 
theory as expressed in terms of orbital true longitude. 


DEFJNITIONS FROM THE MODIFIED HANSEN THEORY 

When the differential equations given in Musen's articles 2, 3 are solved, expressions for the 
components H , T, and of the w function, for the x parameters, and for the perturbation of the pseudo- 
time n 0 Sz result. The functions S,T, and >P are expressed in terms of orbital true longitude and are 
related to osculating elements by the formulas 



T 



e cos <p - 
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e sin 




Here is the deviation of the osculating true anomaly from the true anomaly of the auxiliary ellipse, 
e is the osculating eccentricity, and h is proportional to the reciprocal of the Delaunay variable g; 
that is, 


G 


JL 

h • 


( 2 ) 


The quantities h 0 and e 0 are the elements of Hansen’s auxiliary ellipse. 
The X parameters are defined by the formulas 


X, - sin ~n cos N 


sin ~2 sin N 


a 3 - cos ~2 s ^ n K 


k A = 


cos *2 cos K. 


r (3) 
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Here i is the osculating angle of inclination of the orbit plane and corresponds to I in Brouwer's 
development. The quantities K and N are Fourier series of the order of the perturbations and do not 
contain secular terms. 

The angular variables are given by the formulas 


f - cv - 7 r 0 ~ 4> , 

CD = (g ~ c)v + (tt 0 - 0 Q ) + 4> + K + N 


> ( 4 ) 


6 = (l-h')v + (9 0 +K~N . 


The quantities f , a,, and 6 are the osculating true anomaly, argument of perigee, and longitude of the 
node, respectively. The quantities g, c, and h' in the right hand side of Equations 4 are proportional 
to the mean motions of the argument of latitude, mean anomaly, and the longitude of the ascending node, 
respectively. The quantities w 0 and 8 0 are prescribed constants. 

The time element of the auxiliary ellipse is denoted by the symbol z and is often called the 
pseudo-time. When orbital true longitude is the argument, the mean anomaly of the auxiliary ellipse 
is c(n 0 ) H z . The symbol n 0 appears with different meanings in the articles of Brouwer and Musen; 
therefore the symbol (n 0 ) H is adopted here instead of the n 0 appearing in Musen's article. The quan- 
tity 8z is the deviation of the pseudo-time from the unperturbed satellite time. 


OSCULATING ELEMENTS FOR THE MODIFIED HANSEN THEORY 

By inverting Equations 1 and 3 it is readily found that 

G " h ” h 0 ^ 3 + 27 3 + • • ' ) ’ 

1 /.. 1 (. 

e - e. 


1 if \ 

+ 2 (T“ e 0 H) + 24 (4e 0 S 2 - 4HT j+ . 

/ g H 2 \ 

si - cos i 0 ^l +T + 27 + • • 7 • 

Similarly, the quantities associated with the angular variables are found to be 
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It is instructive to derive an expression for the perturbation of the radius vector of the satellite 
using Equations 5 and 6. This leads to an important result already given in the modified Hansen 
theory. If u denotes the reciprocal of the radius vector of the satellite, and s the deviation of an 
osculating element from its value in the auxiliary ellipse, then 


But we have 


Su = u “ u 


h 0 IE s(£) 
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where f - f - <j>\ and from Equations 5 and 6 we have, to order j 2 , 




8e = 


a 
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T" e 0 H 


Sf = ~ 4 > = 
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> 0 ) 


Substituting the required quantities from Equations 8 into Equation 7 we get, to order J 2 , 

1 h o - H - 

Su = -2 — W +-6“ u , 


( 10 ) 


where w = E + T cos f + f sin f , which is consistent with the results of the modified Hansen theory. 


COMPARISON OF RESULTS TO THE FIRST ORDER IN J 2 

By solving the equations given by Musen, first order analytic solutions for the quantities H , T, 
f, and the X. parameters were obtained by Bailie and Fisher 8 . When the analytic expressions for g, 
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T , and 'P are substituted into Equations 5 immediate agreement is obtained with the periodic part of 
the elements G, e, and I obtained in Brouwer’s solution. Similarly, agreement for the periodic part 
of the expressions for the angular variables & and 6 given by Equation 4 with the variables g and h 
can be readily obtained, when the analytic results of Bailie and Fisher are introduced. 

It has now been indicated that the periodic part of the solution of the elements of the satellite 
problem by Brouwer and Musen agree to the first order in J 2 . Although differences in the arbitrary 
constants and arguments of the trigonometric terms do exist, they do not appear in the first order 
solutions for the trigonometric parts of the elements since they have J 2 as a multiplier. These dif- 
ferences are exhibited in the terms of the second order and are discussed below. 


THE ARBITRARY CONSTANTS OF THE THEORIES 

Differences of order J 2 appear in the arbitrary constants of the solutions of the satellite prob- 
lem by Brouwer and by Musen. The quantities c 0 and c x cos f in Musen’s work are added to the w 
function and consequently the constants c 0 and c % are added to S andT; These constants thus occur 
in the solution for those elements derived from S and T. Constants also appear in the solution for 
the elements by Brouwer. In order to compare the two theories the constants of the elements G, e , 
and cos i with respect to true anomaly are found from Brouwer’s development to order J 2 . Similarly, 
constants with respect to orbital true longitude are found in the article of Bailie and Fisher 8 and 
are listed in Table 1. 


Table 1 


Constants Appearing in the Satellite Theories (Order J 2 .) 


Quantity 

Brouwer's Notation 

Musen's Notation 

(von Zeipel method) 

(modified Hansen method) 
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The constants appearing in Table 1 are defined as follows: 
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= Vl - e " 2 
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These values are taken from the article of Brouwer and from the article of Bailie and Fisher 8 . 

The relations between the mean motions of the argument of perigee and the longitude of the node 
in the articles of Brouwer and Musen are given by the formulas 


ug — f \ 

~St ~ ( n o) H (8 " c > • 

dh w 

“3T = ( n o) H C 1 " h ') • 




> ( 12 ) 




Formulas to order j 2 2 for these mean motions are given in the article of Brouwer and the article of 
Bailie and Fisher. At first sight the terms in j 2 2 seem to disagree. However, by taking the rela- 
tionships given in Table 1 into account, full agreement is obtained to order J 2 2 in the mean motion of 
the variables as defined in Equations 12. 

The differences in the constants given in Table 1 will also be exhibited in the coefficients of 
trigonometric terms of order J 2 2 in the elements derived by the methods of Brouwer and of Musen. 
Additional differences appear in these coefficients and are due to differences in the arguments of the 
trigonometric terms<> We shall now describe these. 


THE TIME ELEMENTS OF THE THEORIES 

In the von Zeipel method adopted by Brouwer the true anomalies f and f ' appear; Brouwer then 
shows how to relate these true anomalies to the true time of the satellite. In the Hansen method 
modified by Musen the true anomaly of the auxiliary ellipse T (or £ as it is denoted by Bailie and 
Fisher 8 ) appears; Musen shows how to relate T to the true time. The true anomalies in the two 
theories differ by trigonometric terms of order J 2 ; consequently, it is logical to apply Taylor’s 
theorem to find the relation between these two true anomalies. 
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We recall that the true anomaly is a function of the eccentricity as well as of the mean anomaly, 
as is shown by the equation of the center. 6 Also, we have 

f = T - <t> 

osc ^ 

e osc = e 0 +Sl 

where f osc and e osc are the osculating true anomaly and osculating eccentricity, whiles and Se are 
of order j 2 . 

The quantity f appearing in Brouwer’s article is related to the osculating mean anomaly i by 
the equation 

37 = i 1 ~ e o 2 ) 2 l 1 + e o cos O 2 • ( 14 > 

We then find by Taylor's theorem that 

F(f) = F(f) + (tf + -|j 8e)|= , (15) 

where 



Bf 

3e 


/ 2 + e Q cos f \ 

\ 1 " e o 2 / 


sin f 


(16) 


Equation 15 transforms a function of the true anomaly of the auxiliary ellipse,?, to a function 
of the true anomaly f appearing in Brouwer's article. To extend this transformation so that a func- 
tion of T may be expressed in terms of the mean true anomaly f ', we simply apply Taylor's Theorem 
again to obtain 

/ Bf df \ 3F . „ 

f<» c • O’) 


to the first order in j 2 . Here i' is the mean true anomaly in the sense given in Brouwer's article 
and may be evaluated by Kepler's equation for a given instant of time. The perturbation A i is the 
deviation of the mean anomaly from its mean value and is given by the formula A/ = - BS^/BL' in 
Brouwer's article. It may also be found from the variation equation in terms of orbital true longi- 
tude by the methods adopted in the article of Bailie and Fisher. 


In particular, if F(f) = sin f , we have 


sin f - sin f ' + \4> + ~ 1 $e + A/) cos f 

V Be (M 


The multiplier of cos f is of order J 2 , so that when P is given f may be found by successive 
approximations. 


7 



It is possible to arrive at Equation 17 by approaching the transformation from a somewhat dif- 
ferent point of view. From the work of E. Brown 9 we have 

„ _ df Bf 

8f " 6 l 8/ + Be 8e ■ (19) 

Since the symbol 8 refers to the deviation of the osculating element from the corresponding ele- 
ment of the auxiliary ellipse, we have 

s l=l- [ c ( n 0 ) H z + /o'] * (20) 

where the quantity in brackets is the mean anomaly of the auxiliary ellipse. Since by Table 1 

c ( n o) h = ‘3T (21) 

and 


z - t Sz , 


Sz 


we have, from Equation 20, 

Sl = j + /„*)]- c(n 0 ) H 

When only the short period terms of the mean anomaly are considered, we have 

61 * tl - c(n 0 ) H Sz . 

Or, substituting in Equation 19, we find 

df . . _ Bf df 

T c K)h 5z - 4 >+ li ?e+ dT A/ > 


( 22 ) 


(23) 


(24) 


which equals the multiplier of BF/Bf in Equation 17. Consequently, Equation 17 may be thought of as 
transforming a function of the pseudo-time z to a function of time t (to order J 2 ) by the relation 


BF 


F(z) - F(t) + ^ Sz . 


(25) 


In order to complete the transformation, the argument of perigee is taken into consideration. 
From Equations 4 it is seen that a term proportional to the equation of the center must be included 
in the transformations. It then follows that 


/ Bf df \ BF 

F(f,*0 = F(f\g # ) +(* + 3 ;- Se A/J + (g - c)(f 


BF 


( 26 ) 


where 


(g - c)v + (w 0 - e 0 ) . 
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If F(f ' ) and F(f ' ,g' ) in the right-hand sides of Equations 17 and 26 represent periodic terms of 
order j 2 , then the subsequent terms in the right-hand sides of these equations will be of order J 2 2 . 
Thus, by means of Equations 17 and 26, it is possible to derive comparisons of the periodic terms 
of the two theories when they are both developed to order j 2 2 . 

Kozai 10 has extended the work of Brouwer 1 to order J 2 2 in the periodic terms and J 2 3 in the 
secular terms. Unfortunately, a corresponding extension has not been made of the analytic results 
of Bailie and Fisher 8 . However, for the purpose of checking the formulas of the present paper the 
author has obtained the development of only the short period terms of the Delaunay variable G to 
order J 2 2 in terms of orbital true longitude. This was done so that comparison could be made with 
the corresponding terms obtained by Kozai using the method of von Zeipel. For this variable, 
agreement was obtained. 


SUMMARY AND CONCLUSIONS 

The solutions to the problem of the near earth satellite without drag given by Brouwer 1 and by 
Musen 2 - 3 agree when carried out to the same order in J 2 . Due allowance must be made for the dif- 
ferences in the constants and in the ways of expressing the time element. 

The differences in the arbitrary constants have been tabulated here to the first order in j 2 , 
and transformations have been given relating the true anomaly of the auxiliary ellipse to that of 
the satellite. 


REFERENCES 

1. Brouwer, D., "Solution of the Problem of Artificial Satellite Theory Without Drag," Astronom. J. 
64(9):378-397, 1959. 

2. Musen, P., "Application of Hansen’s Theory to the Motion of an Artificial Satellite in the Gravi- 
tational Field of the Earth," J. Geophys. Res. 64(12):2271-2279, December 1959. 

3. Musen, P., "The Theory of Artificial Satellites in Terms of the Orbital True Longitude," 

J. Geophys. Res. 66(2):403-409, February 1961. 

4. Kozai, Y., "The Motion of a Close Earth Satellite," Astronom. J. 64(9):367-377, November 1959. 

5. Garfinkel, B., "The Orbit of a Satellite of an Oblate Planet," Astronom. J. 64(9):353-367, 
November 1959. 

6. Brouwer, D., and Clemence, G., "Methods of Celestial Mechanics,” New York: Academic Press, 
1961. 

7. Bailie, A., and Bryant, R., "Osculating Elements Derived from the Motion of an Artificial Satel- 
lite," Astronom. J. 65(8):451-456, October 1960. 


9 



8. Bailie, A., and Fisher, D., "An Analytic Representation of Musen’s Theory of Artificial Satellites 
in Terms of the Orbital True Longitude," NASA Technical Note D-1468, January 1963. 

9. Brown, E., "An Introductory Treatise on the Lunar Theory," Oxford: Cambridge University 
Press, 1896. 

10. Kozai, Y., "Second-Order Solution of Artificial Satellite Theory Without Air Drag Astronom. J. 
67(7):446-461, September 1962. 


10 



Appendix A 


Analytic Results from the Modified Hansen Theory 


By solving the equations given by Musen*, first order analytic solutions for the quantities «, T, 
ip } } \ 3 f gj i -h', and g -c were obtained by Bailie and Fisher* . These are listed here for the con- 

venience of the reader. 

To conform to modern notation, the quantities k 2 , k 3 , and k 4 appearing in the modified Hansen 
theory are here designated byj 2 /2, j,/2, and -J 4 /8 respectively. In formulas A1 through A5 the 
quantities g and c introduced by integration have been restored, along with the quantity m ; these were 
set equal to unity by Bailie and Fisher who were then concerned only with terms of order J 2 . 


S = J h ^ (1-3 cos 2 i 0 ) 2 - ^ - e 0 J 


9 h o 4 / „ x J cos 2r) c os (j + 27?) , _ cos - 2??) 

- T J 2 ~^2 l 1 “ cos x o) \ i~ + e 0 C + 2g + e 0 2g - C 
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10 e 0 2 cos 2 i 0 
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. 3 J3 h Q 2 . 

(2g - 27,) J. e o sin 1 o Sln ^ 


15 il V 

•* £ T C 


16 J 2 ^ c o 


8 cos 4 if 


1-3 cos 2 i 


0 1-5 cos 2 i r 


cos (2£ - 2rj) ; 


(Al) 


1 ho 4 


?= 8 J 2— 2 (1 - 3 cos2 L o) 


+ 6e 0 cos 2g + e Q 2 cos 3£ 


8e„ + 4e 0 3 

1 Ja V , 2 \ 
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1 +Jl - e 0 2 J 

8 c fj? r 13 c oj 


^J 2 i 1 ~ cos2 i o) { Se 


cos 2rj 
0 g 


♦Musen, P., ’’Appl ication of Hansen’s Theory to the Motion of an Artificial Satellite in the Gravitational Field of the Earth,” J. Geopbys. 

tBlilie^A^and 7 Fish« ,’ D.^^AnAnaTytic Representation of Musen’s Theory of Artificial Satellites in Terms of the Orbital True Longitude, 
NASA Technical Note D- 1468 , January 1963 . 
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